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First-principles calculations are applied to study spin Hall effect in semiconductors and simple 
metals. We found that intrinsic spin Hall conductivity (ISHC) in realistic materials shows rich sign 
changes, which may be used to distinguish the effect from the extrinsic one. The calculated ISHC in 
n-doped GaAs can be well compared with experiment, and it differs from the sign obtained from the 
extrinsic effect. On the other hand, the ISHC in W and Au, which shows opposite sign respectively, 
is robust and not sensitive to the disorder. 
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The existence of Berry phase in systems with spin- 
orbit coupling (SOC) can act as gauge field in the mo- 
mentum space, which in turn affects the transport behav- 
ior of electrons in real space, and produces the fascinat- 
ing new phenomena in solid crystals. One typical exam- 
ple is the intrinsic spin-Hall effect (ISHE) proposed re- 
cently 0,0- The ISHE is an effect that, for non-magnetic 
materials with SOC, a transverse pure spin transport can 
be induced by external electric-field in the absence of 
magnetic field (even at room temperature). It is dis- 
tinguished from the extrinsic spin-Hall effect (ESHE), 
which is due to impurity scattering |3(. The obvious ad- 
vantages of ISHE, especially for the field of spintronics, 
have stimulated extensive studies recently, both theoret- 
icaUy |iiSMp^IIlIIlIIlIIllllIIlII3 and 
experimentally [lj|, Ujj ■ Up to now, most of the theo- 
retical studies were done based on certain model Hamil- 
tonians, like Luttinger, Rashba, or Dresselhauss Hamil- 
tonians. These studies provide pictures for the under- 
standing of deep physics (such as the effects of vertex 
correction 20]), but on the other hand, neglect the band 
details, which could be very important (as will be ad- 
dressed in the present paper) due to the topological na- 
ture of ISHE. Two experimental evidences have been pro- 
vided for the existence of spin Hall effect (SHE) HlHil. 
The SHE on the 2D hole gas HI] is likely of the intrin- 
sic origin, however, the intrinsic or the extrinsic origin 
of the SHE in the 3D electron film ^| is still under de- 
bate 0, H3| ■ In order to have close comparison between 
theory and experiment, parameter- free considerations in- 
cluding all band details are highly desirable, and will be 
the main focus of this paper. In particular, our calcula- 
tions make reliable predictions on the sign of the ISHE, 
which in some cases differs from the sign obtained from 
the ESHE. This qualitative difference can be used to de- 
termine the origin of the effect. 

In this letter, we will consider realistic materials by us- 



ing first-principles calculations to study the ISHE in var- 
ious systems, including semiconductors (Si and GaAs), 
and simple metals (W and Au). The main difficulty of 
such study comes from the accurate evaluation of Berry 
curvature. We have recently developed a technique to 
evaluate such property accurately, and have applied it 
to the calculation of the anomalous Hall effect (AHE) 
in Fe and SrRuOs [23I E3 | . Here we apply it to study 
the ISHE. Besides the quantitative evaluation, we will 
concentrate on the sensitivity of ISHE to band details 
and the rich sign changes of ISHE in various materials, 
which will provide strong support for future experiments 
to identify the ISHE. 

First-principles calculations have been done based 
on standard density functional theory using accurate 
FLAPW (full potential linearized augmented plane wave) 
method, in which the relativistic SOC has been treated 
fully self-consistently. The exchange-correlation poten- 
tial was treated by the generalized gradient approxima- 
tion (GGA), whose validity for the systems we consider 
here has been shown by many other studies. The under- 
estimated band gap in GGA for semiconductors does not 
produce problems for our purpose here, because the con- 
cerns are for the d.c. limit (w=0). Accurate k-point in- 
tegration has been done by tetrahedron method or adap- 
tive mesh refinement [13 ■ The convergence of calculated 
results with respect to the number of k-point has been 
carefully checked. In general, the number of k-point re- 
quired to achieve accuracy of 5% is around 1,000,000 in 
the Brillouin Zones (BZ). For all the calculations pre- 
sented here, experimental lattice parameters are used. 

Suppose the external electric field is applied along the y 
direction, then the linear response of the spin (a z ) current 
along the x direction can be obtained from the Kubo 
formula by evaluating the spin Hall conductivity tensor, 
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where \nk > is the eigen wave function of Bloch state 
with eigen value E n and Fermi occupation number 
v y is the velocity operator, j x is the spin current operator, 
which is defined as j(o~ z v x + v x a z ). fl„ is the spin Berry 
phase connection of the Bloch state, and is responsible 
for the anomalous transverse transport we studied. The 
important point here is that all band details and SOC 
are self-consistently taken into account (no adjustable 
parameters), no approximation beyond linear response 
theory has been used. It is also straight forward to take 
into account the impurity scattering effect by allowing 
finite life time broadening S, and the finite temperature 
effect in the Fermi distribution f n f.. For the definition of 
sign, positive a xy means that spin-up (s z =l/2) compo- 
nent flows to the positive x direction. For the convenience 
of comparison, in our following discussions, we convert 
spin conductivity into the unit of charge conductivity by 
multiplying a factor ^ip to the calculated values. 
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FIG. 1: The calculated (a) band structure; (b) spin-Hall 
conductivity a xy as function of Fermi level position for bulk 
GaAs. We define the converged Fermi level without doping as 
energy zero point, and evaluate a xy by rigidly shift the Fermi 
level position. The panel (c) gives the a xy of n-GaAs as func- 
tions of carrier (electron) density after subtracting the part 
that does not contribute to spin accumulation (see the text 
part for explanation). Note the factor of 10 for the T=30K 
and 5=16meV curve. 

Semiconductors: Fig.l and Fig. 2 show the calculated 
band structures and a xy as functions of Fermi level (Ep) 
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FIG. 2: The calculated (a) band structure; (b) spin-Hall con- 
ductivity a xy as function of Fermi level position for bulk Si. 



position for GaAs and Si respectively. In the following, 
we will first concentrate on the zero-temperature and 
clean limit (T=0K, <5=0eV), then the effects of disorder 
and elevated temperature will be addressed later. 

For hole-doped GaAs (Fig. 1(b)), the calculated a xy is 
large, it reaches about 300 ri" 1 cm _1 , which is about the 
same order of magnitude as that estimated from Lut- 
tinger model 0. For hole-doped Si (Fig. 2(b)), on the 
other hand, the maximum of a xy is about 50 f2 -1 cm -1 . 
This can be understood by taking into account the fact 
that the strength of SOC in Si is about 1/7.7 of that in 
GaAs. In both cases, the o~ xy first increases with increas- 
ing hole density, after reaches maximum it goes down. By 
carefully calculate the Berry curvature distribution il(k), 
we found that the contributions from the light hole and 
split-off bands are mostly negative, which will compen- 
sate the positive contributions from the heavy hole band, 
and finally suppress the o~ xy for large hole density. 

For the electron-doped GaAs (Fig. 1(b)), the situation 
is complicated and far beyond what we understood from 
model analysis. The calculated a xy shows sign changes as 
electron density varies: negative (positive) for low (high) 
density. Such behavior is related to a small splitting of 
the conduction band due to the lack of inversion symme- 
try in GaAs. For n-GaAs, the conduction band bottom 
has mostly s orbital character. Due to the s-p hybridiza- 
tion, some p characters exist in this band, leading to the 
Dresselhaus type SOC in the simplified model (therefore 
the vertex correction is not important [20lE3 |l. 

For electron-doped Si (Fig. 2(b)), on the other hand, 
the obtained o~ xy at clean limit is quite large (although 
the strength of SOC in Si is small), and has negative 
sign. In this case, the conduction-band-bottom is nei- 
ther around T point nor s orbital like. The calculated 
Berry phase contributions to ISHE are mostly related to 
the conduction bands around X point of the BZ. Consid- 
ering the factor that the spin relaxation time in electron- 
doped Si is typically much longer than that in hole-doped 
case |23, our results suggest the possibility to realize 
ISHE in Si, which is the most important semiconductor 
material. 
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For both GaAs and Si, the effects of temperature 
and disorder are important for the electron-doped cases, 
while not so dramatic for hole-doped ones. As shown 
in Fig.2(b) for Si, by putting <5=20meV and T=300K in 
our calculations, the a xy is significantly reduced for n-Si, 
while not so much for p-Si. This is also true for GaAs 
as shown in Fig. 1(c). The disorder and temperature will 
even cause sign change for a xy in n-GaAs as will be dis- 
cussed later in comparison with experimental results. 

For the insulating GaAs and Si (Ep located in the 
gap), however, the calculated a xy is non- vanishing (about 
43 and 7 ft -1 cm -1 for GaAs and Si respectively). First 
we have to emphasize that this is not due to numerical 
error, which is four orders of magnitude smaller. Our re- 
sults can be regarded as a, generalization of the concept 
of "Spin Hall Insulator" Murakami et.al. studied 

the spin Hall effect in narrow-gap and zero-gap semi- 
conductors like PbTe and HgTe, which have "special" 
band structures (such as the inverted light hole and heavy 
hole bands in HgTe), and demonstrate the existence of 
spin Hall insulator. However, the results here suggest 
that such "special" band structure is not necessary in 
general. In real materials, there always exist finite hy- 
bridizations, which will produce non- vanishing ISHE in 
insulators. Note that ISHE is not quantized [lj, in qual- 
itative difference with the AHE. Nevertheless, we should 
emphasize that the existence of such ISHE in insulator 
will not produce any spin accumulation due to the lack 
of broken time reversal symmetry [l6| . 

To make comparison with experimental results on n- 
doped GaAs [lj|, we show in Fig. 1(c) the calculated a xy 
as functions of electron density, with the subtraction of 
the part that does not contribute to the spin accumu- 
lation (the value within the gap). Now it is very clear 
that the a xy (for T=0K, S=0meV) is negative for small 
doping, but change sign to be positive for large dop- 
ing. We also notice that such fluctuation of a xy are 
suppressed by introducing disorder and temperature ef- 
fects. This, on the one hand, is the nature results of 
topological origin of ISHE, and on the other hand, sug- 
gests complication of ISHE in realistic materials. For 
the experimental doping density (3 x 10 16 cto~ 3 ) 0], the 
calculated a xy (<5=0meV, T=0K) has the same sign (neg- 
ative) as that obtained in experiment, and also agrees 
with Ref. This is in sharp contrast with the ex- 

trinsic SHE, which has opposite sign as estimated in 
Ref. The absolute value of calculated a xy at clean 
limit is two order of magnitude larger than experimental 
value (-0.005ft _1 cm _1 ), however a compatible number 



issue remains to be checked in the future. 



(-0.01ft" 



_1 ) can be obtained by introducing a finite 



lifetime broadening S=2meV as shown in Fig. 1(c). Unfor- 
tunately experimental parameter 6 for unstrained sample 
is unclear 19] . Using the measured p xx (300 ft/mi) for 
strained sample, we estimate the <5=16meV, which give 
positive a xy (T=30K, see Fig. 1(c)). Nevertheless, con- 
sidering the uncertainty of experimental parameters, this 
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FIG. 3: The same notation as shown in Fig. 2, but for W. 
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FIG. 4: The same notation as shown in Fig. 2, but for Au. 

Simple Metal: The ISHE in simple metal has not been 
studied yet, although it is not surprising to expect that 
ISHE exists in such systems, due to the same mechanism. 
We chose elemental W and Au as examples because of 
the relatively larger SOC. For W the charge conductivity 
mostly comes from 5d states around Fermi level, while 
for Au it is mostly from 6s states. The calculated band 
structures and a xy are shown in Fig. 3 and Fig. 4 for W 
and Au respectively. Besides the very strong ISHE ob- 
tained (the a xy can reach as high as -1390 ft _1 cm _1 for 
W, and 731 ft _1 cm _1 for Au), we notice that a xy is neg- 
ative in W, but positive in Au. This again suggests the 
rich sign changes of ISHE in realistic materials. What is 
more interesting is that the ISHE in W and Au are robust 
and not sensitive to the disorder (in opposite to GaAs or 
Si). By adding a very large disorder effect (5=0. 5eV), 
the calculated a xy only change slightly. Given these spe- 
cial characters, we suggest that both W and Au are nice 
candidates for future experimental examination of ISHE. 
Especially for Au, where the conduction electrons have 
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mostly s characters, relatively long spin relaxation can 
be expected |25j 

Sign Issue and Discussions: As presented in our above 
results, the ISHE shows very rich sign changes, which 
are independent of the carrier type and the sign of im- 
purity potential: (1) for Si, the sign of ISHE is the same 
as ordinary Hall effect, i.e., positive for hole doping and 
negative for electron doping; (2) for W and Au, however, 
the sign of ISHE is opposite with their carrier type: W 
(Au) has hole (electron) type conductivity but negative 
(positive) ISHE; (3) for n-doped GaAs, the sign of ISHE 
changes with increasing doping. Such rich sign changes 
are more than what we can expect from extrinsic scat- 
tering mechanism Two mechanisms, namely skew 
scatterin g an d side-jump, were mainly discussed in lit- 
eratures [23. For a simple discussion, we consider the 
skew scattering mechanism ,27], which dominates over 
side-jump contribution for the weak disorder limit. In 
this case, the si gn o f ESHE depends on the sign of scat- 
tering potential 27] . It is natural to expect that the sign 
of skew scattering does not change with changing impu- 
rity density. However, we predict that the sign of ISHE 
can change with the same type of doping (in n-GaAs). 
Such difference may be used to distinguish ESHE from 
intrinsic contributions. We can also use simple metal W 
and Au for such purpose, because opposite signs of ISHE 
(which has nothing to do with scattering potential) are 
predicted. Nevertheless, the sign issue should be regarded 
as a very important aspect of SHE, and can be used in 
future experiments. 

In summary, we have performed detailed studies on 
the ISHE for various realistic materials based on accu- 
rate parameter-free first-principles calculations, and pre- 
dict rich sign changes of ISHE. Furthermore, we demon- 
strate that the ISHE in semiconductors (GaAs and Si) is 
highly sensitive to band-details and disorder, while ISHE 
in simple metals (W and Au) is robust and not sensitive 
to disorder. The calculated ISHE for n-GaAs can be well 
compared with experimental results, while the extrinsic 
spin Hall contribution has opposite sign. 
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